Dissipative continuous Euler flows in two and three dimensions 



A. Choffrut, C. De Lellis, L. Szekelyhidi Jr. 
March 3, 2013 

Abstract 

In |DLS12 two of these authors construct dissipative continuous (weak) solutions 
to the incompressible Euler equations on the three-dimensional torus T 3 . The building 
blocks in their proof are Beltrami flows, which are inherently three-dimensional. The 
purpose of this note is to show that the techniques can nevertheless be adapted to the 
two-dimensional case. 



Contents 

1 Introduction 

2 Setup and general considerations on the construction 4 

2.1 Construction of (v\,pi, R±) 6 

2.2 A linear set of stationary flows 

2.3 The geometric lemma 

3 The maps v\ , p\ , R\ |ni 

3.1 The perturbation w Q 1C 

3.2 The constants r] and M 12 

3.3 The correction w c 13 

3.4 The pressure term p\ 13 

3.5 The tensor R\ [13 

4 Estimates 

4.1 Schauder estimates 15 

4.2 Estimates on the corrector and the energy 16 

4.3 Estimates on the Reynolds stress [2C 



1 



5 Conclusions |23 

5.1 Proof of Proposition [3] 23 

5.2 A quantified version of Proposition [3] 24 

5.3 Proof of Theorem M 



1 Introduction 

In this paper we will take 

d = 2 or 3 

and we will consider the incompressible Euler equations 

dtv + div (v <g) v) + Vp = 0, div u = . (1) 

The following Theorem was proved for d = 3 in the paper [DLS12] and we will show in 
this note that: 

• the same statement holds for d = 2 and indeed the proof of [DLS12] can be suitably 
modified to handle this case as well; 

• in d = 3 the approach yields nonetheless solutions which are genuinely 3-dimensional 
(cp. with Theorem [2j). 

Theorem 1 Let T > 0. Assume e: [0, T] — > R is a smooth positive function. Then, there is 
a continuous vector field v. T d x [0, T] — > M. d and a continuous scalar field p: T d x [0, T] — > R 
solving the incompressible Euler equations in the sense of distributions and such that 

e(t)= \v\ 2 (x,t)dx, t£[0,T\. (2) 

If (v,p) is a C 1 solution of (HJ, we can scalar multiply the first equation by v and use 
the chain rule to derive the identity 

a^ + div,((^ +P ),)=o. 

Integrating this last equality in space we then derive the conservation of the total kinetic 
energy 

— / \v\ 2 (x,t)dx = 0. (3) 
dt Jjd 

Thus, classical solutions of the incompressible Euler equations are energy conservative and 
are therefore "rigid" compared to the continuous solutions, for which Theorem Q] shows that 
any energy profile e(t) is indeed possible. 
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The existence of continuous solutions which violate the conservation of the total ki- 
netic energy was first suggested in |Ons49j by Onsager, where indeed he conjectured the 
existence of such solutions in 3 space dimensions with any Holder exponent smaller than 
|. Onsager also asserted that such solutions do not exist if we impose the Holder conti- 
nuity with exponent larger than | and this part of his conjecture was proved in |Eyi94| 
and |CET94j . The considerations of Onsager are motivated by the Kolmogorov theory of 
isotropic 3-dimensional turbulence, where the phenomenon of anomalous dissipation in the 
Navier-Stokes equations is postulated. This assumption seems to be widely confirmed ex- 
perimentally, whereas no such phenomenon is observed in 2 dimensions. Indeed, for d = 2 
the conservation law for the enstrophy does prevent it for solutions which start from suf- 
ficiently smooth initial data. However, the considerations put forward by Onsager which 
pertain the mathematical structure of the equations do not depend on the dimension and 
this independence appears clearly also in the proof of [CET94] . which works for any d > 2. 

The first proof of the existence of a weak solution violating the energy conservation 
was given in the groundbreaking work of Scheffer [Sch93| . which showed the existence of a 
compactly supported nontrivial weak solution in R 2 x R. A different construction of the 
existence of a compactly supported nontrivial weak solution in T 2 x R was then given by 
Shnirelman in |Shn97j . In both cases the solutions are only square summable as a function 
of both space and time variables. The first proof of the existence of a solution for which 
the total kinetic energy is a monotone decreasing function has been given by Shnirelman in 
[ShnOOj . Shnirelman's example is only in the energy space L°°([0, oo[, L 2 (R 3 )). 

In the works [DLS09[ IDLS10] these existence results were extended to solutions with 
bounded velocity and pressure and the same methods were also used to give quite severe 
counterexamples to the uniqueness of admissible solutions, both for incompressible and 
compressible Euler. Further developments in fluid dynamics inspired by these works ap- 
peared subsequently in |Chil21[UFGlH[ShvllllSzell|lSWlHIWiellj and are surveyed in the 
note [DLS11 . The paper [DLS09] introduced a new point of view in the subject, highlight- 
ing connections to other counterintuitive solutions of (mainly geometric) systems of partial 
differential equations: in geometry these solutions are, according to Gromov, instances of 
the /i-principle. It was also observed in [CDLSll] that the Onsager's Conjecture bears 
striking similarities with the rigidity and flexibility properties of isometric embeddings of 
Riemannian manifolds, pioneered by the celebrated work of Nash |Nas54| . It turns out that 
the iteration procedure leading to Theorem [T] shares some fundamental aspects with the 
approach of [Nas54| : we refer to the introduction of |DLS12| for a thorough discussion. 

As already mentioned, we will also show that the 3D flows constructed in Theorem Q] 
are genuinely three-dimensional. In order to formulate our statement precisely, consider 
a solution (v,p) of ([I]) on T 3 x [0, T] and denote with the same letters the corresponding 
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solution on M 3 x [0, T] with the obvious periodicity in space. The solution is then not 
genuinely 3-dimensional if, after suitably changing coordinates in space, it takes the form 



v(x,t) = (v 1 (x 1 ,X2,t),v 2 (x 1 ,X2,t),v 3 ) (4) 

where V3 is a constant. Now, a careful analysis of the proof of Theorem [1] will give the 
following statement. 

Theorem 2 Let d = 3. For any e > the solutions constructed in Theorem [TJ can be 
chosen to satisfy 

< £, (5) 



sup 

te[o,i] 



v(x, t) ® v(x, t) dx e(t)ld 

T 3 3 



and 



SU P IMIif-i(T 3 ) < e. (6) 
*e[o,i] 



It is then easy to conclude that, for sufficiently small e, a solution as in Theorem [2] cannot 
have the form (jlj). Indeed, from the bound (|5|) we would have |t?3| > gy/min tg [ 0j i] e(t), 
whereas from dSJ) we would have \vs\ < e. 

Remark A similar analysis in the two-dimensional case leads to the conclusion that the 
two-dimensional flows constructed in Theorem Q] are typically not parallel flows. However, 
it is classical that such flows are necessarily stationary, and this is not possible if e(t) is not 
constant. 



2 Setup and general considerations on the construction 

As already mentioned, the proof of Theorem Q] is based on an iteration procedure in each 
step of which one solves a system of equations closely related to the Euler equations. We 
introduce some terminology. We let S dxd denote the space of symmetric d x d matrices 
and SQ Xd its subspace of trace-free matrices. Assume v,p,R are smooth functions on T d x 
[0, T] taking values in R d ,M, and SQ Xd respectively. We say that (v,p,R) solves the Euler- 
Reynolds system if it satisfies 

dtv + div (y ® v) + Vp = div R, div v = 0. (7) 

Theorem [1] is a consequence of the following 

Proposition 3 Let e: [0, T] — > R be a smooth positive function. Then, there exist constants 
r] and M , depending on e, with the following property. 

Let 6 < 1 and (v,p,R) be a solution of the Euler- Reynolds system 

dtv + div (v <8> v) + Vp = div R 
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such that 



3 r 5 

-8e(t)<e(t)- \v\ 2 (x,t)dx < -6e(t), te[0,T] (8) 

4 J T d 4 

and 

sup\R(x,t)\ < r)8. (9) 

Then there exists (v\,pi, R\) solving the Euler- Reynolds system and satisfying the following 
estimates: 

3 f 5 

-8e(t) < e(t) - / \v l \ 2 (x,t)dx < -6e(t), t e [0,T] (10) 

8 y 8 

sup \Ri(x,t) < -7)8, (11) 
x,t * 

sup |«i (x,t)-v(x,t)\<My/S, (12) 

and 

sup|pi(x,t) -p(x,t)\ < MS. (13) 

Proof of Theorem [TJ Start with vo = 0, po = 0, i?o = 0, and 8 = 1. Apply 
Proposition [3] iteratively to obtain sequences v n , p n , and R n solving the Euler- Reynolds 
system (J7J) and satisfying 



^—<e(t)-[ \v n \ 2 (x,t)dx < t 6 [0,21, 

4 2" " u hd y J "4 2" L J 

sup|i? n (x,t)| < — 

x,t 2" 



sup|t> n+ i(a;,t) - v n (x,t)\ < M \J^, 

. , . .. M 
Sup \p n+ x{x,t) -p n (x,t)\ < — . 
a;,t 2 n 

The sequences t> n and p n are Cauchy in C(T x [0, T]) and hence converge (uniformly) 
to continuous functions v and p respectively. Likewise, R n converges (uniformly) to 0. 
Moreover, taking limits in the estimates on the energy, 

\v\ 2 {x,t)dx = e(t), te[0,T]. 

Also, we may pass to the limit in the (weak formulation of the) Euler-Reynolds system ([7]) 
and this shows that v,p satisfy the (weak formulation of the) Euler equations (pQ). ■ 

In the rest of this Section we collect the main ingredients in the proofs of Theorems Q] 
and H 
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2.1 Construction of (vi,pi,Ri) 

Assume (v,p, R) as in Proposition [3l Write v\ = v + w and p\ = p + q where w and q are 
to be determined. Then, 

dfVi + div (v\ ® vi) + Vpi = div (tt; ® to + qld + i?) + c^w + div (o^ioltogo). 

The perturbation w should be so chosen as to eliminate the first term in the right-hand 
side, viewing the remainder as a small error. Note how this first term is reminiscent of the 
stationary Euler equations. Roughly speaking, the perturbation w will be chosen as a high 
oscillation modulation of a stationary solution. More specifically, it will be taken of the 
form 

w = w + w c 

where w Q is a highly oscillator term given quite explicitly and the corrector term w c en- 
force the divergence-free condition divu; = 0. The oscillation term w Q will depend on two 
parameters A and p such that 

A,m,- e N. 

V- 

In fact, A and p will have to be chosen sufficiently large, depending on v, R, and e, in order 
that the desired estimates (fTUj) . (fTTj) . (fT2j) . and (fT3|) be satisfied. 



2.2 A linear set of stationary flows 

The stationary Euler equation is nonlinear. The following Proposition says that there exists 
a linear space of stationary solutions. We first introduce some notation. 

• Case d = 2 For k G Z 2 , we let 

so that 

&fc(0 = V^fc(e), where = (-%,%) 
hence div^ = 0. Observe also that frj. = and ip^ = ip-k- 

• Case d = 3 For /c G Z 3 , we let 

&*(£) := B fc e lfc -« (15) 

where will be chosen appropriately, see Lemma In particular, one should have 
div^ bk = and = 
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Lemma 4 (A linear set of stationary flows in 2D) Let v > 1. For k G I? such that 
\k\ 2 = v, let dfc € C suc/i £/iat a& = a_fc. Then, 

w(o = e a kh(o, *co = E a ^(£) 
ifei 2 =i/ ifci 2 =v 



are W-valued and satisfy 



div^W ® = V 5 ( Li- - ^ ) . (16) 



Furthermore, 

In other words, W and \I/ are the velocity field and stream function, respectively, of a 
stationary flow with pressure P = — — h i^-^-. 

Proof By direct computation one finds A^ipk = — \k\ 2 ipk, and hence that A^ = — z/\E'. 
Recall the identities 

div^W ® W) = \W\ 2 + (curl^W) W ± 

where curl^W = d e W 2 - d^W 1 = A^ and W 1 - = (-W 2 , W 1 ). Then, 

\W\ 2 

as desired. 

As for the average, write 

f W®W(OdC = ^ajakbj^bk 
Jt 2 77 

bl 1*1 



E ajOfce 



If j / k, f T2 e*(*-'H = and it is 1 if j = k. Thus, 

ln,J 2 I TH - 

1*1 1*1. 

in-r— " 

where the last identity follows from prj" ® prj" = fid — jlj ® by direct calculation. 
The following three-dimensional version was proved in [DLS12] . 



|fc| 2 =l> 11 11 |fc|2=J/ V 
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Lemma 5 (Beltrami flows) Let v > 1. There exist Bj, G C 3 /or |fc| = v such that, for 
any choice of at G C such that~a~k = a-k, the vector field 

W(0 = a kh(0 (19) 

\k\=v 

is divergence- free and satisfies 

\W\ 2 

div (W®W) = V- — -. (20) 



Furthermore, 

£ W®Wd£ = 

|*|=i/ 
the 

Proof This is Proposition 3.1 in [DLS12] . 



^^^(id-A^A). (21) 



Note that in the three-dimensional case the pressure is given by P = — ^S-. 



2.3 The geometric lemma 

Lemma 6 (Geometric Lemma) For every N G N we can choose ro > and v >\ with 
the following property. There exist pairwise disjoint subsets 

f {keZ 2 : \k\ 2 = u} if d = 2 
AjC \ { * G Z3 : |^=, } if d = 3 ' 

and smooth positive functions 

7? G C°°(i? ro (Id)), j G {1, . . . , JV}, G A, 

swc/i £/iai 

i. A; G Aj implies — k G Aj and 7^ = 7^2; 
£. /or eac/i G -B ro (Id) we /iaue i/ie identity 

R = E U' ^)) 2 ( Id - m ® ) R e ^o(W). (22) 

Proof This was proved for the case d = 3 in Lemma 3.2 of |DLS12j (up to a normalization 

constant). Here we treat the case d = 2 only. For each »6l 2 \ {0} we define 

v v 
M v = Id - — <g> — . 

\v\ \v\ 

Step 1 Fix u > 1 and for each set F C {A; G Z 2 : |A;| 2 = z^} denote c(-F) the interior of the 
convex hull in 5 2x2 (the space of symmetric 2x2 matrices) of the set Mp := {M& : k G F}. 
We claim in this step that it sufices to find v and N disjoint subsets Fj C {k G 1? : |A;| 2 = u} 
such that 
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• c(Fj) contains a positive multiple of the identity. 

Indeed, we will show below that, if Fj satisfies these two conditions, then we can find ro > 0, 
a subset Fj C Fj, and smooth positive functions a[ G C°° (i?2r (Id) ) for k G Tj such that 

We then define the sets Aj and the functions 7^ for k G Aj according to 

• \0) = if k g Aj \r j; 

•7? == J^U^? forked 

Note that the sets Aj are then symmetric and that the functions 7^ satisfy ()22|) . Moreover, 
since at least one of A^l is positive on B2 ro (ld), 7^ is smooth (and positive) in B ro (ld). 

We now come to the existence of Tj. The open set c(Fj) contains an element aid for 
some a > 0. Observe that the space 5 2x2 of symmetric 2x2 matrices has dimension 3. Since 
aid sits inside the open set c(Fj), there exists a 4-simplex S with vertices A±, . . . , A4 G c(Fj) 
such that aid belongs to the interior of 5. Let then $ so that the ball U centered at aid 
and radius 1? si contained in S. Then, each R G U can be written in a unique way as a 
convex combination of Ai's: 

4 

R = J2p i (R)A i 

i=l 

where the functions are positive and smooth in U. 

Using now Caratheodory's Theorem, each A{ is the convex combination ^ n \ n M v . n 
of at most 4 matrices M Vin where Uj >n G Fj, where we require that each Aj >n be positive. 
(Caratheodory's Theorem guarantees the existence of 4 elements M Vi n such that A{ belongs 
to their closed convex hull. If we insist that all coefficients should be positive, then some of 
these elements should be thrown away.) 

Set now ro := Then 

R = Y -Pi(aR)X i>n M Vin , R G B 2ro (Id) 
L — ' a 

i,n 

and each coefficient ^/3j(ai?)Aj jn is positive for R G £2^0 (Id). The set Tj is then taken as 
{vi tn }. Since one might have Vi >n = v\^ m for distinct (i,n) and (l,m), the function A& will 
be defined by 

X k (R) = V -Pi(aR)X ijn 
L — ' a 
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and this completes Step 1. 

Step 2 By Step 1, in order to prove the lemma, it suffices to find v and N disjoint families 
• • • , Fjy C y / i/S 1 nZ 2 such that each set c(Fj) contains a positive multiple of the identity. 
Note that S 1 n Q 2 is dense in S 1 . Indeed, let 

™ t \ / 2it u 2 — 1\ , m n 
u£ R^ s (u) ~ J E S C M 2 . 

yu^ + 1 u z + 1 J 

Clearly s(Q) C Q 2 . Since Q is dense in R and s is a diffeomorphism onto S 1 \ {(0, 1)}, the 
claim is proved. 

In turn, there exists a sequence Uf. — > +oo such that the sets S 1 n "^=^ 2 converge, in 
the Hausdorff sense, to the entire circle S . Given the sequence v^, one can easily partition 
each y^S 1 n Z 2 into N disjoint symmetric families {iylj^i,...,^ m such a way that, for 
each fixed j, the corresponding sequence of sets |^=i^| converges in the Hausdorff sense 
to S 1 . Hence, any point of c(S 1 ) is contained in c{^=Fj?) for k sufficiently large. On the 
other hand, it is easy to see that c(§ 1 ) contains a multiple aid of the identity. (One can 
adapt for instance the argument of Lemma 4.2 in |DLS09| .) 

By Step 1, this concludes the proof. ■ 



3 The maps 

Let e(t) be as in Theorem [1] and suppose (v,p,R) satisfies the hypothesis of Proposition [3l 
In this Section we define the next iterates v% = v + w, pi = p + q, and R\. The perturbation 
w will be defined as w = w + w c where w is given by an explicit formula, see (|30l) , and 
the corrector w c guarantees that div w = 0. 

3.1 The perturbation w Q 

We apply Lemma [6] with 

N = 2 d 

(where d E {2, 3} denotes the number of space dimensions) to obtain v>\ and ro > 0, and 
pairwise disjoint families A» with corresponding functions j¥ E C°°(B ro (Id)). 

The following is proved in |DLS12j and works for any number of space dimensions. 
Denote C\,... ,Cn the equivalence classes of Z d / ~ where k ~ I if k — I E (2Z) rf . The 
parameter /j£N will be fixed later. 
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Proposition 7 (Partition of the space of velocities) There exists a partition of the 
space of velocities, namely M>-valued functions cei(v) fori € 7L d satisfying 

£(a,(«)) a = l (23) 



such that, setting 

4'W) = E a <(He _i(fe >, j = l,-..,N, k G Z d , (24) 
we have the following estimates: 

S up\D™$\v,T)\ < C(m,d)n m , (25) 

V,T 

sup iDHSir^ +*(*;■ t)| < CKIfcM)^ 1 . (26) 



(7) (7) 
Furthermore, 4>^ = (p_ k and 



l*?W)l a = £"?(«)■ (27) 



Set now 

^-wri'-rl 1 " 121 *'^) (28) 

and 

fl(j/,*):=p(a)Id-fl(i/,a) (29) 
(observe that the tensor — R is then the traceless part of R). Define 

w Q (x, t) := W(x, t, Xt, Xx) (30) 

where 

W(y,s,T,£) = ^2 a k(y,s,T)b k (C) 
|fc| 2 =i/ 

= v%jEE^(^)#w^).*K)' ( 31 ) 



j=i fceA 



(The velocity fields bk were defined in (114p for d = 2 and in (|15j) for d = 3.) For d = 2 we 
also introduce the corresponding stream function 

ip (x, t) := t, At, Ax) (32) 
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where 



\k\ 2 =u 

(The stream functions ipk were defined in (|14p .) 



3.2 The constants 77 and M 

In this Section we fix the values of the constants r\ and M from Proposition [3l 

The perturbation w Q is well defined provided — G i? ro (Id) where ro is as in Lemma [BJ 
By definition ()28[) of /? and assumption ([8]), 



> ~r; — w~re(i) > c<5 min e(t) 



where c= M^F' Then ' 



Id 



P(*) 

Thus, we choose r] satisfying 



< 



R 



< 



c5 min 4e[0jT] e(i) ' cmim, e[0)T] e(i) ' 



n < -c min e(t)ro- (33) 
2 te[o,T] 



Observe that this restriction is independent of 5. 
We choose first a constant M' > 1 such that 



M ' > 4 E E sW^)! supl^Cv.r)] sup|6 fe (C)| (34) 

V^lfcGA, R 
and then choose M > 1 such that 



,sup n< . <T e(t) 
4d(2vr) d 

Observe that with these choices we have 



■s/M 

max{v / ^||V'o||o, H^ollo} < — 7T~ V\J. (35) 



since < ^yr<5e(t) by (EHJ) and 
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3.3 The correction w c 

To obtain w from w we need to introduce the Leray projection onto divergence-free vector 
fields with zero average. 

Definition 8 (The Leray projector) Let v G C°°(T d ,M. d ) be a smooth vector field. Let 



Q v ■= + f v 

where (p G C°°(T d ) is the solution to 

A(f> = div v in T d subject to 

We denote by 



V:=L-Q 

the Leray projector onto divergence- free vector fields with zero average. 
The iterate v± is then expressed as 

v\ = v + Vw = v + w + w c , w c = —Qw = w — w . (36) 



3.4 The pressure term pi 

We set 

' V- (^-"f) if d = 2, 
Pi-=< I l 2 J 37 

p-W- if d = 3. 

This choice for p\ will become clearer at the end of the proof, see Lemma [T51 

3.5 The tensor R x 

To construct R\, we introduce another operator. 

Definition 9 Let v G C°° (T d ,M. d ) be a smooth vector field. We define IZv to be the matrix- 
valued periodic function 

ryu + (v«r-(div«)id if d = 2 

{ \ (VVu + (VVu) 1 ) + § (Vu + (Vu)- 1 ) - |(div u)Id if d = 3 
where u G C 00 ^,^) is tfie solution to 

Au = v - f v mT d , subject to -f u = 0. 

The operator 7?. satisfies the following properties. 
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Lemma 10 (K = div^ 1 ) For any v £ C°°(T d ,R d ) we have 

1. TZv(x) is a symmetric trace-free matrix for each x 6 T d ; 

2. div TZv = v — / Td v . 

Proof The case al = 3 is treated in |DLS12j and thus we assume d = 2. Clearly, TZv is 
symmetric. Next, 

tr TZv = 2div u — 2div u = 



and 



div TZv = An + div (Vu) T - div ((div n)Id) 
= An + div(Vn) T - Vdivn 
= An 



v — -f v 
It 2 

by definition of n. ■ 
Then we set 

i?i := TZ (d t vx + div (vt <g> v{) + Vpi) . (39) 

One verifies, as in |DLS12| (after Lemma 4.3, p. 15), that the argument in the right-hand side 
has zero average: div (v\ ®V\ +pild) clearly has average zero, and so does dtv\ = dtv + dtw 
since dtv = —div [v v + pld) has average zero as well as w by definition of V . In turn, 
Lemma [10] yields 

dtV\ + div (v± ® vi) + Vpi = div R\. 

4 Estimates 

The letter m will denote a natural number (in N), and a a real number in the interval (0, 1). 
The letter C will always denote a generic constant which may depend on e,v,R,i/,a, and 
S, but not on A nor n. We will further impose that 

1 < n < A. (40) 

The sup- norm is denoted ||/||o = sup T d |/|. The Holder seminorms are given by 

[f) m ■= max H^/llo, 
\j\=m 

[f\m+a ■■= max sup — — — 

h\=m Xf Ly \x — y\ 
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and the Holder norms are given by 



j=0 



m+a 



+ [/], 



We also recall the following elementariy identity: for < r < 1 

[fg] r <C([f] r \\g\\o + \\fh[g] r ). 

4.1 Schauder estimates 

The next Proposition collects estimates in Holder spaces ( "Schauder estimates" ) for various 
operators used in the remainder. 

Proposition 11 For any a G (0,1) and m G N i/iere exists a constant C = C(m,a,d) 
satisfying the following properties. If 4>,ip: T d — > R are the unique solutions to 



Acf> = f 







= divF 



then 



\m+2,a <c||/ll m,a. 
Moreover, we have the following estimates: 



and 



i+l, a 

<C\\F\ 



\\Qv\\ m +a 
W'P'V 1 1 m+a 
||^||m+l,a 
||^div A|| m+Q 
|7£Qdiv A|| m+Q 



^ C|M|m+a, 

— C 1 1 v 1 1 m+a j 

— C 1 1 v 1 1 m+a j 

— C||^4||m+a, 

— C||^4||m+a- 



(41) 
(42) 
(43) 
(44) 
(45) 



These estimates are fairly standard and a detailed proof is given in [DLS12] in three dimen- 
sions. It is easy easy to see that they also hold in two dimensions as well: the only operator 
which is defined differently is 1Z and this makes the 2-dimensional case only easier. Suffice 
it to say that these estimates are the expected ones: V and Q are differential operators 
of degree 0; 1Z is a differential operator of degree —1; and div is a differential operator of 
degree 1. 

The effect of the oscillation parameter A is described in the following 



Proposition 12 Let k G Z d \ and A > 1 be fixed. 
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1. For any a G C°°(T d ) and m G N we have 

[ a(x)e iXk - x alx 

2. Let 4> x G C°°{T d ) be the solution to 

A0 A = / A in T 



- \m ■ 



d 

subject to f T d4>\ = where f\(x) = a{x)e lXk ' x — j Td a{y)e tXk ' y dy . Then, for any 
a G (0, 1) we have the estimate 

||V0 A |U < ^iz^ll«llo + + -^[a] m+a 

where C = C(m,a,d). 

This was established in [DLS12] and is in fact valid in any dimension. 

The following is a consequence of the definition (|38|) of 1Z, the Schauder estimate (f4~2]) 
for V, and Proposition 1121 

Corollary 13 (Estimates for the operator 1Z) Let k G Z d \ be fixed. For a smooth 
vector field a G C°°(T d ,R d ), let F(x) := a{x)e iXk ' x . Then, we have 

c c c 

for C = C (m, a, d) . 

4.2 Estimates on the corrector and the energy 

We recall that w Q (x,t) is defined in (130|) . In what follows, if / is a function of time and 
space, we will use the notation ||/|| r for the "Holder norm in space", that is 

11/11,= SUp \\f(;t)\\r- 

te[o,T] 

Lemma 14 1. Let a k G C°°(T d x [0,T] x R) be as in (EZJ). Then, for any r > 0, 



\\ a k(- 


s. 


r)\\r 


< 




\\d s a k (- 


s. 


r)\\r 


< 


Cfi r+1 


\\d T a k {- 


s. 


r)\\r 


< 




\\{d T a k + i(k ■ v)a k )(- 


s. 


r)\\r 


< 


cy -1 
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2. The matrix- valued function W ® W is given by 

(W®W)(y,s,T,Z) = R(y,s)+ U k (y, s, r)e lk < (46) 

l<\k\<2v 

where the coefficients U k G C°°(T d x [0,T] x S dxd ) satisfy, for any r > 0, 

\\U k (;s,r)\\ r < C//, 
\\d s U k (;s,r)\\ r < CfJ, r+1 , 

\\d T U k (;S,T)\\ r < C(jT, 
\\(d T +i(k-v)U k )(-,S,T)\\ r < Cff- 1 . 



3. For d = 2, we have 

^-4- E ^ 

l<|fe|<2i/ 

where the functions d k satisfy for any r > 

||ofc(-,s,r)|| r < C//'. 

In a// i/iese estimates the constant C depends on r, e, v, and R, but is independent of 
(s, r) and fj,. 

Proof The first two items are proved in |DLS12j for d = 3 and we only briefly recall their 
proof since it is identical for d = 2. The estimates on the coefficients a k follow from the 
estimates (|25j) and (j26|) on <p^ ■ Next, write W <S) W as a Fourier series in £: 



W®W(y, S ,r,e) = ^o(y,s,r)+ ^ U k (y,s,T)e ik - 

K\k\<2u 



where the entries in the U k s are quadratic in the a k s and ||afc||o < C. Thus, the L^'s satisfy 
the claimed estimates. Furthermore, 



U (y,s,T) = i W®Wd£ 



j=i fceAj zeC, 

i?. (47) 
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As for the third item of the Proposition, concerning — u%- when d = 2, we compute, 



2 "2 

2 _ IZJ2 



omitting variables and remembering that the sums are over j and k such that \ j\ = \k\ = u, 

j,k 



^E a i a * 



V 



j,k 



i{k-j)-t 



j,k 

v ^-^ v — ' 

\k\*=V J^fc 

J E h 2 + E ^^) gl(fe - 

|fc|2 = jy l<|fc|<2i/ 







where the coefficients a fc are quadratic in the a k s. But from the expression (|17p for (T4 / (8>W / )^ 
we deduce 

|T^| 2 =tr(W®W) = tri?+ ^ tr[/ fc e ifc '« 

l<|fc|<2i^ 

= E h 2 + E trC/ ^ K - 

|fc| 2 =v l<|fc|<2i/ 

Subtracting the above expression for u^ 2 from that of |Vt^| 2 we obtain the desired expression 
with suitable coefficients a k which satisfy the same estimates as the a k s. ■ 



Lemma 15 (Estimate on the corrector) 

IKIla < (48) 

Proof This is proved in Lemma 6.2, p. 21, in [DLS12j for d = 3. For clarity we reprove it 
for d = 2 with the appropriate adjustments. Recall that 

w Q (x,t) = V] a k (x,t,\t)V ± i) k (\x) = Vf$>(x,t,\t,\x) 

\k\ 2 =u 

where *(y,s,r,£) = Y^\k\ 2 =u a k{v, s , r)ip k {0- But 

V 1 (a fc (x,t,At)^ fe (Ax)) = Aa fc (x,t,At)(V ± ^)(Aa;) + ^(Ax)V ± ci fc (a;,t,At) 



or 



a fe (x, i, At) (V -1 ^) (Az) = ^{ v± («fc(a;,i,At)^(Ax)) - Vfc(Ax)V ± a fc (x,t,A)} 
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and thus 



w (x,t) 



J V± Yl a k (x,t,\t)MXx)\ - \Y1 ^fc(Ax)V ± a fe (x,t,At). 

\|jfc|2=i/ J \k\ 2 =u 



Since Q eliminates the divergence-free part and div o V = 0, we have by (|36|) 

w c (x,t) = Qw (x,t) = J] ^ fc (Aa;)V- L afc(a;,t,At) =: -Qu c (x,t). (49) 

Thus, by Schauder estimate (|4T|) for Q and the estimates on the coefficients a k from 
Lemma [Til we find 



C 



(I 



— \ 11 \l-a 



A 1 



Lemma 16 (Estimate on the energy) 

e(t)(l- 6 -)- [ \ Vl \ 2 dx 

2 JT d 



< c- 



A 



l-a ' 



(50) 



Proof This is proved in Lemma 6.3, p. 21 of |DLS12j for d = 3. For clarity we briefly 
recall it in the case d = 2. Taking the trace in the expression (|46p for TV <g> W gives 

|W| 2 = tr J R+ ^ trU k e ik <. 

l<\k\<2v 

From part (1) of Proposition [12] with m = 1, and using estimates on from Lemma [1 
we find 

This, along with the estimate (|48|) on u> c and ||to ||o — C, implies 



w Q \ 2 — tr R dx 


and 


/ v ■ w dx 


< 











I 1 2 I 1 2 I |2 j 

t>i — u — \w \ ax 



< c- 



l-a ' 



Now by definition (|28l) of p we have tr R = dp = ,,A a (e(t) (l — | ) — J Td |t;| 2 dx). Putting 



the above together finishes the proof. 
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4.3 Estimates on the Reynolds stress 



In order to clarify the choice for p\ as in (|37p . we will temporarily write p\ = p + q. With 
this, we have 

div Ri = dfVi + div (vi i>i) + Vp\ 
= d t w + v ■ Vw 

+div (w (3>w + qld + R) 

+dtiv c + div (v\ ® w c + w c §<) v\ — w c ® w c + v ® w Q ). 

We split the Reynolds stress tensor into the transport part, the oscillation part, and the error 
as shown on the right-hand side of the above identity. In the remainder of this Section we 
estimate these terms separately. 



Lemma 17 (The transport part) 



li + A 1 



\\K(d t w + vVw )\\ a < C[ — + ) . (51) 



Proof This is proved in Lemma 7.1, p. 22 of [DLS12j for d = 3 and the proof given there 
is valid for d = 2 as well. ■ 



Lemma 18 (The oscillation term) 



2 

(div (w ® w + qld + R)^j <Cti^- ( 52 ) 



n 

Proof This is proved in Lemma 7.2, p. 23 of |DLS12j for d = 3. The main difference in 
the case d = 2 is in the role of q. 
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Recalling that R(y, s) = p(s)ld — R(y, s), see (f29]h and noting that p = p(t) is a function 
of t only, 

div (w ® w + R + qld) 
= div (w <3 w — R + qld) 

= div(. ^. -^-(M-4)id) + v(, + M-4 

\w\ 2 ^ 2 

d\vy[W ®W - R - [ J-J- - v— ) Id 



\W\ 2 ^ 2 
Adiv f I W ® W - I J—i- - i^— I Id 



2 2 

^2 
— f 

2 

/ , K| 2 <Ao 

/ /|w|2 ^ 2 

div, f TV (gi W - R - f - i/— ) 

div, f ]T C/ fc e iAfc -- £ S fe e iAfe ^Id 

\|fc| 2 =i; \k\ 2 =u 

where two cancelations occur by construction of w , see (1161) . and by definition of q, see 
(|37p . The estimate follows from Corollary [13] with m = 1 . ■ 



Lemma 19 (The error - I) 



\\Kd t w c \\ a < C^. (53) 



Proof This is proved in Lemma 7.3, p. 23 of |DLS12j for d = 3. For clarity we briefly 
recall it for d = 2 with the appropriate adjustments. 

Recall that u c = ^|fc| 2 =^ ' l l ; k(^x)V ± ak{x, t, Xt) was defined in (f39j) and thus 

d t u c (x,t) = A ^ k {Xx)V L d T a k {x,t,Xt) 

\k\ 2 =u 

+ ^2 V , fc(Ax)V- L (? s a fc (x,t, At) 

|fc| 2 = !/ 

But for any vector- valued function A(y,s,r), we have 

A(x, t, At) ® j^e^J = »AA(x, t, ^)^p + e lAfc ' x ( ^3 • Vj A(z, i, At) 

hence 

M^x)A(x, t, \t)e lXk - x = ^div (a(x, t, Xt) ^3 e iAfc ") - ^ Xk ' x (J^ ■ v) A(x, t, Xt). 
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Therefore dtu c is of the form 

d t u c = div U c + u c 

where 

\\U c \\ a < CfiX a , \\u c \\ a < Cfi 2 X a 
owing to the estimates from Lemma [TH In turn, 

\\Kd t w c \\ a < j(\\TZQdivU c \\ a + \\TZQu c 
C 

< -r (\\U c \\ a + \\u c \\ a ) 

< 

A 

owing to the estimates ([H]), (@3|), and P5|) . 



Lemma 20 (The error - II) 

||ft(div {vi ®w c + w c ®v 1 -w c ® w e ))\\ a < Ctt^- ( 54 ) 

Proof This is proved in Lemma 7.4, p. 24 of [DLS12] for d = 3 and the proof given there 
is valid in the case d = 2 as well. ■ 



Lemma 21 (The error - III) 

||fc(div (v <g> w ) 



< C 



A 1- 



(55) 



Proof This is proved in Lemma 7.4, p. 24 of [DLS12j for d = 3 and we briefly indicate 
the proof in the case d = 2. Using that div^ = 0, we find 

div (v w ) = w ■ Vv + (div w )v 

= ^2 a k (bk -V)v + (Va k ■ b k )v 



|fc| 2 =i/ 

E 

\k\ 2 =u 



<~lk 



(ik 1 



- 



Vl*l 



• V }v + Va k ■ 



ik 1 - 

w 



J,\k-x 



The estimate follows from Corollary 1131 with m = 1. 
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5 Conclusions 



5.1 Proof of Proposition [3] 

Recall that e(t) is given as in Proposition [3] and that (v,p,R) is assumed to solve the 
Euler- Reynolds system (J7J) and to satisfy the bounds © and Q. 

We have now all estimates available in order to fix the parameters fi, A, and a so that 
the estimates (|10p . (jlip . (|12|) . and (|13p may hold. For simplicity we will take 



/' 



\? (56) 



for some (3 to be determined (although strictly speaking this can only hold up to some 
constant depending only on (5 since it is required that /x G N). Recall that C denotes a 
generic constant (possibly) depending on e, v, R, v, a, and 8, but not on A nor \i. 
Recall that the constant M of Proposition [3] has already been fixed in (|35p so that 



max{^||^ ||o, |K,||o} < — ^— v 7 ^. (57) 



2 

Since v\ — v = w + w c , the bound (|12[) on di - « follows provided 

IKL < C-£- = CA"^ 1 < V5. 
~ A — 

The bound (fT3|) on pi — p = — ( z/% ) also follows from ([57|) . 



2 " 2 ^ 

The bound (|10p on the energy follows from (|50p provided 

C-^- = CA^" 1 < - min e(t). 
A 1 -" ~ 8 te[o,T] K ' 

Finally, the estimates $TQ, §2$), (JSHJ, flM}, and (J55J, imply that 

||Ai||a < C(A a_/3 + A^" 1 + A 2 ^" 1 ). 
In conclusion, imposing 

a < /3 and a + 2/5 < 1 (58) 



ensures that the bounds (|T0|) . (|TTj) . ([12]) . and ([13]) hold provided A is chosen sufficiently 
large. 

This concludes the proof of Proposition [3l 
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5.2 A quantified version of Proposition [3] 



In what follows we fix d = 3. The proof of Theorem [2] follows from a closer analysis of the 
proof of Theorem [I] in particular from a more precise version of Proposition [3] which will 
be stated below. In order to prove it we collect first a series of estimates which make some 
of the statements in the proof of Proposition [3] more precise. 

Lemma 22 Let e,v,p,R be given. Let R,w ,w c ,vi,Ri be defined as in (QSfl), $EW> 
and A39\). Then, there exist constants C = C(e,v,p, R) and 7 > such that 



\\W c \\a 

v ■ w Q dx 



T 3 



w <8> w Q dx — R dx 



v 1 <8> v\ — v ® v — w <S> w dx 



la 



< CA" 7 

< C\~~ f 

< C\~~ f 

< CA" 7 . 



(59) 
(60) 

(61) 

(62) 
(63) 



\\R 

Proof Take < a < (3 < 1 satisfying a + 2(3 < 1, see ([35]) . and set 

7 = min (1 — a — 2/3, /3 — a, 1 — 2a — (3) . 

Estimate (j59[) is Lemma H3J Estimates (I60p . (I6ip and (I62D follow from the proof of 
Lemma HH Specifically, (160p follows since u is fixed and u; is oscillatory, and (1611) fol- 
lows since w ® w — R \s the sum of oscillatory terms. For (|62j) we write 



v 1 (8> «i dx 



IP 



T 3 



v®vdx+ w ®w dx+ I (v ® w + w ® v + w ® w c + w c ® w ) dx. 



The estimate follows since zu is oscillatory and it? c is small (u is fixed) . Estimate (|63p follows 
from Lemmas 7.1-7.4 of [DLS12] , ■ 



Proposition 23 Let e be as above and let v,p,R solve the Euler- Reynolds system ([7]). 
Suppose that there exist < 5 < 1 satisfying 



e(t)(l-5) 



\v(x, t)\ 2 dx 



< Z e® 



\\R\\o < rjS. 

Let R, vi = v + w = v + w + w c , p\ = p + q, R\ be defined as in §EEj) , Wfy > 

and \39\) . so that they satisfy the Euler- Reynolds system: 

d t v\ + div (v\ g) Vi) + Vpi = div R\. 
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For any e' , the following inequalities are satisfied provided A is chosen sufficiently large 
depending on (v,p,R). 



*>HH. 



\vi(t)\ 2 dx 



\Vl - v\\ 



SUP \\Vl - V\\ H -1( T 3) 

t 



< oe(*) 



v\ <8> vi dx 



v vdx 



w ® w dx 

'T 3 

w ® w — Rdx 



r 3 



< 


min ( 


< 


MV6, 


< 


e', 


< 


M5 


< 


e' 


< 


e'. 



Proof With 



we have 



P(t) 



1 



3(2vr) 3 
min <t<r e(t)8 



<t) (i - - f M*>*)\ 2 dx 

m&x <t< T e(t)6 



12(2tt) 3 

Owing to the definition (|33|) of 77 we obtain 



< P(t) < 



4(2vr) J 



(64) 

(65) 

(66) 
(67) 

(68) 

(69) 

(70) 

(71) 
(72) 




















p 



< r . 

Thus, vi,pi,Ri etc. can be defined and estimated as in Lemma [22j We will now choose A 
in Lemma [221 sufficiently large (depending on e,v,p,R) so that the desired estimates hold. 

Recalling the definition (j29[) of R in terms of p, we have 



e(t) [1-"-) - J \vi(t)\ 2 dx 
e(t) (l-0 - f^\v{t)\ 2 dx ■ 
- / (\w \ 2 -trR) dx+ (\vi 

JJ3 Jj3 



\w \ dx + 



T 3 
2 



(Kl 



w, 



I'p 
\w a \ 2 ) dx. 



o| 2 ) dx 



Invoking (|6ip and (|62p . we obtain (|64p with large A. 

From (|63p it is clear that the bound (|65p on i?i holds with A large. 

From (|59p . w c can be made arbitrarily small with A large, while from (|30p and (|34[) we 
have |tt> | < ^M'^/p < \M\fb. Thus, the bounds (i66j) on v\ — v = w G + w c and f[68|) on 
Pi — p = — 2 follow with A large. Furthermore, w; is a fast oscillating function and w c 
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is small in the sup-norm, so that we achieve the H 1 -bound (|67p on v\ — v as well with A 
large. 

Finally, flBU) and fl62]) trivially imply ([69]) and (170]) with large A. ■ 



5.3 Proof of Theorem [5] 

We are now ready to prove that the solution v,p constructed in Theorem Q] in the case d = 3 
satisfies the estimates of Theorem [2] provided A is chosen sufficiently large at each iteration. 
Set 

1 ■ ,s 
£n := - mm em. 

u 2 o<t<T y ' 

By choosing s > smaller, if necessary, we may assume that e < |eo- AVe recall that the 
solution (v,p) is obtained as a limit 

v:=limv^ n \ p : = limp^ 

n n 

where the sequences v^ n \ p^ are as follows. We let 

v m = o, p <°) = o, #°) = o. 

For n > 0, we construct 

v (n+l) = v (n) + ^(n+1) ? p (n+l) = p (n) + g (n+1) > £(n+l) 

using Proposition [23l with = v, p^ = p, R^ n+1 ^ = R\, t?( n+1 ) = v\, p( n+1 ) = p 1 

5 = 2~ n , and e' = e2~ n - 1 , (73) 

so that 

d t v {n) + div («W ® + Vp (n) = divi? (n) . 

Recall that, in Proposition 1231 the velocity v\ is obtained by adding two perturbation w 
and w c to u, defined, respectively, in (|3"U]) and in (J3SD - We use the notation w^ 1 ^ = w Q 
and w;i n+1 ^ = w c . Moreover, we set 

,W(t) = ^ (e(t) (1 - 2—1) - 

and i?( n ) = p( n )ld — R^ (observe that therefore p^ n > = p and R^ = R, with p and as 
in (J2SJ) and ([29"]) ). In particular, 
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according to ([71]) with v = and 5=1 (step n = 0). 
Now 

oo 

|H| H -i<Elk (n+1) l|H-i<e 



n=0 



easily follows from (|67p . This establishes 
Next, for each n, 



T 3 

+ / 

/T3 



+ 



+ 
+ 



(v {n+1) ® ^ {n+1) - w (n) ® « (n) - ^ n+1) ® wi n+1) ) da; 
+ (2vr) 3 p (n) Id - / dx 
u^ n+1) ®^ n+1) -i? (n) ) dx 



v (n+l) <g w (n+l) _ ^(n) ^ v (n) _ ^(n+1) ^ ^(n+lA ^ 



But the sum of the last two integrals is bounded by e2 n using ([M]) . ([70]) . and ([75]) . There- 
fore, taking the trace we obtain 



v( n+1 \x,t)\ 2 dx = / \v^ n) {x,t)\ 2 dx + ^fp^ n \t) + 0{e2^ n ) 



T 3 



e(t)(l-2- n - 1 ) + 0(e2- n ), 



where 0(e2 n ) denotes an error term bounded by e2 n . In particular this yields p^ n \t) 
_L_^2~ n_1 e(i) + 0(e2~ n ). Furthermore, by summing over n we obtain 



/ 



v <S> v dx 



oo „ oo 

/ R [n) dx + (2vr) 3 P (n)ld + °( £ )- 

n=0'' Ti n=0 



Since, by ([65]) and ([73]) we also have 



oo „ 

y / RWdx 



< e 



we deduce 



Jt 



v($vdx = -e(i)Id + 0(e) . 
3 3 
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